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Abstract 

Background: Both fear and pain processing are altered in post-traumatic stress disorder (PTSD), as evidenced by 
functional neuroimaging studies showing increased amygdala responses to threats, and increased insula, putamen 
and caudate activity in response to heat pain. Using psychophysiology and functional magnetic resonance 
imaging, we studied conditioned and unconditioned autonomic and neuronal responses in subjects with PTSD 
versus trauma-exposed non-PTSD control (TENC) subjects. A design using an electric shock selected by subjects to 
be 'highly annoying but not painful' as an unconditioned stimulus (US) with partially reinforced cues allowed us to 
partly disentangle the expectancy- and prediction-error components from sensory components of the 
unconditioned response. 

Results: Whereas responses to the conditioned stimulus (CS) were similar in PTSD and TENC, the former displayed 
higher putamen, insula, caudate and amygdala responses to the US. Reactivity to the US in the anterior insula 
correlated with PTSD symptom severity. Functional connectivity analyses using the putamen as a seed region 
indicated that TENC subjects had increased amygdala-putamen connectivity during US delivery; this connection 
was disengaged in PTSD. 

Conclusions: Our results indicate that although neural processing of fear learning in people with PTSD seems to 
be comparable with controls, neural responses to unconditioned aversive stimuli in PTSD seem to be increased. 



Background 

In classic fear conditioning, an initially neutral, and later 
conditioned stimulus (CS+) is paired with an aversive 
unconditioned stimulus (US) such as an electric shock. 
After pairing, the neutral stimulus is able to elicit fear 
on its own. Much has been learned about the neural 
responses induced by the CS+ presentation. In healthy 
people, conditioned stimuli activate the amygdala, brain- 
stem, insula and parts of anterior cingulate cortex 
(ACC) [1]. Less is known about changes occurring in 
response to the presentation of the more biologically 
relevant US [2]. Responses to a US in healthy subjects 
are accompanied by increased activity in the brainstem, 
thalamus, and cingulate, sensory and insular cortices 
[3-5], structures also known to respond to noxious sti- 
muli [6]. 
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Understanding the neural mechanisms that mediate 
aversive unconditioned responses may be clinically rele- 
vant. Brain regions involved in fear learning, such as the 
amygdala and the ACC, are implicated in the pathophy- 
siology of several anxiety disorders such as post-trau- 
matic stress disorder (PTSD). Other intrinsically 
aversive stimuli such as trauma reminders, aversive 
images and fearful faces also elicit increased responses 
in several brain regions in PTSD, such as the amygdala 
and the ACC [7]. The unconditioned response repre- 
sents a specific exemplar of a more general class of 
aversive stimuli to which people with PTSD may be 
more sensitive. In addition, understanding how aversive 
somatosensory stimuli are processed in PTSD may help 
to explain why pain processing is altered in this disorder 
[8], and clarify the high comorbidity between PTSD and 
chronic pain [9,10]. Recent functional magnetic reso- 
nance imaging (fMRI) studies of pain processing in 
PTSD reported increased insular and putamen reactivity 
to painful heat [11,12]. However, it is unknown whether 
such alterations in pain processing are driven by 
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anticipatory or perceptual properties of noxious stimula- 
tion, and whether aversive but not overtly painful soma- 
tosensory stimuli are also processed differently in PTSD. 

In the present study, we examined the neural correlates 
of responses to electrical shock stimuli in subjects with 
PTSD and in trauma-exposed non-PTSD (TENC) control 
subjects. Within an fMRI fear-conditioning paradigm, we 
examined neuronal and autonomic responses to CS+ and 
US presentations in a partial-reinforcement paradigm 
[13]. The US was a 'highly annoying but not painful' elec- 
tric shock delivered after 62.5% of the CS+ presentations. 
This experimental design allowed us to compare neural 
responses to both the CS+ and US, and to the anticipa- 
tion and unexpected omission of the US (unreinforced 
CS+ trials). By contrasting time intervals representing 
delivered versus non-delivered US, the sensory compo- 
nent of a US response may be in part be isolated from its 
expectancy- and prediction-error related components. 
Moreover, we related alterations in brain responses eli- 
cited by the US to PTSD symptom severity. Lastly, in a 
post hoc analysis, we explored inter-regional functional 
connectivity differences in US-elicited responses in the 
regions in which the subjects with PTSD were found to 
differ significantly from the TENC subjects. 

Methods 

The Partners Healthcare System Human Research Com- 
mittee approved the study, and all subjects gave 
informed consent. 

Subjects 

In total, 23 right-handed people with PTSD and 28 
right-handed healthy TENC subjects were recruited 
from the community. All were administered the Struc- 
tured Clinical Interview for the Diagnostic and Statisti- 
cal Manual, fourth edition (DSM-IV) [14] to determine 
PTSD diagnostic status and possible comorbid axis I 
disorders. 

For the TENC group, participants with any current 
mental disorder were excluded. Subjects with PTSD who 
had current substance dependence were also excluded, as 
were subjects who had used any psychotropic medication 
within 4 weeks before participation (1 year for neurolep- 
tics). Data from four PTSD and four TENC subjects were 
excluded because of excessive head movement while in 
the scanner, leaving a total sample size of 19 PTSD and 
24 TENC subjects. Data from a subset of these partici- 
pants have already been reported in studies examining 
fear extinction [15,16]; however, their responses to the 
US have not previously been reported. 

Rating scales 

The Clinician-Administered PTSD Scale (CAPS) [17], 
Neuroticism-Extroversion-Openness Personality 



Inventory-Revised (NEO PI-R) [18], Spielberger State 
Trait Anxiety Inventory [19], Beck Anxiety Inventory 
(BAI) [20] and Beck Depression Inventory (BDI) [21] 
were administered to all participants. These results, 
along with demographic data, type of trauma, age at 
first trauma exposure, and current comorbid disorders 
are shown in Table 1. 

Fear-conditioning procedure 

We used a partial-reinforcement classic conditioning 
paradigm that has previously been described in detail 
[13]. In the fMRI scanner, each trial began with an 
image of a room (context) containing an unlit lamp, 
which was presented for 3 seconds. The lamp was then 
switched on to show one of three colors (blue, red or 
yellow). Two of the colors (CS+) were followed by an 
electric shock (US) in 62.5% of the cases. The third 
color was never followed by a shock (CS-). The illumi- 
nated lamp was presented 32 times, giving a total of 16 
safe trials (CS-), 10 CS+ trials followed by the shock, 
and 6 CS+ trials in which the shock was omitted. Each 
CS presentation lasted for 6 seconds. Between trials, a 
black screen was displayed for 12 to 18 seconds. The 
lamp color sequence was counterbalanced across sub- 
jects in pseudo-random order (see Additional file 1 Fig- 
ure SI). 
Electric shock 

The US consisted of a 500 ms train of 1 ms spikes at 50 
Hz delivered to the second and third fingers of the right 
hand, with currents ranging from 0.2 to 4.0 mA. Before 
the experiment, the shock current was individually 
adjusted so that it would be perceived as 'highly annoy- 
ing but not painful' by the individual patient. Adminis- 
tration of the US immediately followed CS offset. 
Skin conductance 

Skin-conductance responses (SCRs) were measured 
using a skin-conductance coupler (S71-23; Coulbourn 
Modular Instruments, Allentown, PA, USA) with 8 mm 
(sensor diameter) radiolucent Ag/AgCl electrodes (Bio- 
Pac Systems Inc., Goleta, CA, USA). Electrodes were 
filled with isotonic paste, and placed on the palm of the 
subject's left hand. The skin-conductance electrodes 
were separated by 14 mm. The SCR for each CS trial 
was calculated by subtracting the mean skin-conduc- 
tance level 2 seconds before CS onset from the highest 
skin-conductance level recorded during the following 
CS duration of 6 seconds. SCR responses during the 
interval after the US, the omitted (non-delivered) US or 
the CS- offset were calculated by subtracting the mean 
skin-conductance level recorded during the first 2 sec- 
onds of this interval from the highest skin-conductance 
level during the ensuing 3 seconds. Because autonomic 
responses are delayed with respect to their antecedent 
neural activity [22], this 2-second baseline was more or 
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Table 1 Demographics, comorbidities and type of trauma exposure in the cohort 3 (significant P values in bold). 



PTSD b TENC C P 

Demographics 

Gender, M/F 9/1 0 (n = 1 9) 11/13 (n = 24) 0.5 d 

Age,±SD 36 ± 1 2 (n = 1 9) 30 ± 12 (n = 24) 0.16 

Years of education 1 5 ± 2 (n = 1 9) 1 7 ± 6 (n = 1 9) 0.1 
Personality 6 

NEO neuroticism 27±8(n = 16) 17±9(n = 22) 0.001 

NEO extraversion 22±8(n = 16) 30 ± 6 (n = 22) <0.001 

NEO openness 28±6(n = 16) 32 ± 5 (n = 22) 0.04 

NEO agreeableness 29±6(n = 16) 31 ± 6 (n = 22) 0.3 

NEO conscientiousness 33 ± 7 (n = 16) 34 ± 7 (n = 22) 0.66 
Clinical measures 

CAPS f score 67 ± 24 (n = 19) 8 ± 9 (n = 24) < 0.001 

BDI g score 21 ± 12 (n = 19) 2 ± 3 (n = 24) < 0.001 

BAI h score 1 9 ± 1 3 (n = 1 7) 5 ± 8 (n = 22) < 0.001 

Spielberger Trait anxiety 54 ± 12 (n = 18) 35 ± 9 (n = 21) < 0.001 

Spielberger State anxiety 44 ± 12 (n = 18) 34 ± 12 (n = 22) < 0.001 

Age at trauma, years 24 ± 1 5 (n = 1 9) 1 9 ± 1 2 (n = 1 9) 0.23 

Time since trauma, years 11±15(n = 19) 12 ± 12 (n = 19) 0.95 

Subjects with childhood' trauma n 9 5 0.1 0* d 
Type of trauma exposure, n 

Motor vehicle accidents 3 4 1.0 d 

Sexual assaults 8 3 0.04 d 

Physical assaults 5 7 1.0 d 

Childhood abuse 6 3 0.1 5 d 

Combat 3 0 0.07 d 

Witness to trauma 3 7 0.47 d 
Current comorbidities, n 

Major depression 5 0 0.01 d 

Panic disorder 3 0 0.07 d 

Alcohol abuse 1 0 0.44 d 

Other substance abuse 2 0 0.1 9 d 

Eating disorders 2 0 0.1 9 d 



a Not all ratings were available for all participants; numbers are given in brackets. Participants might also have had more than one type of traumatic event or 
comorbid disorder. 
b Post-traumatic stress disorder. 
c Trauma-exposed normal control. 

d Fisher's exact probability test. All other data were analyzed with Student's f-test, two-tailed. 
e As measured by the Neuroticism-Extroversion-Openness (NEO) Inventory 
f CAPS, Clinician-Administered PTSD Scale. 
g BDI, Beck Depression Inventory. 
h BAI Beck Anxiety Index. 
'Before the age of 16 years. 

less uncontaminated by the response to the shock. A 
representative example of SCR is provided (see Addi- 
tional file 1 Figure S2). 

Image Acquisition 

MRI data were collected using a 3.0 Tesla whole-body 
MRI system (Trio; Siemens Medical Systems, Iselin, NJ, 
USA) equipped for echo planar imaging with a 12-chan- 
nel head coil. Subjects were instructed to lie as still as 
possible, and head movement was restricted with foam 



cushions. After an automated scout image was obtained 
and automated shimming procedures performed, a high- 
resolution, Tl-weighted, three-dimensional, magnetiza- 
tion-prepared rapid acquisition gradient echo 
(MPRAGE) volume was collected to facilitate spatial 
normalization and positioning of the subsequent scans. 
Functional MRI images, sensitive to blood oxygenation 
level-dependent (BOLD) contrast, were acquired with an 
interleaved gradient echo T2*-weighted sequence (TR = 
3000 ms, TE = 30, flip angle = 90°), collected in 45 
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coronal oblique slices tilted 30° down from the anterior- 
posterior commissural line. The voxel size was 3 ± 3 ± 
3 mm. 

fMRI data analysis 
Preprocessing 

SPM8 (Wellcome Trust Center for Neuroimaging, 
http://www.fil.ion.ucl.ac.uk) was used to process all MRI 
data. Structural images were segmented and spatially 
normalized to the Montreal Neurological Institute 
(MNI) 305 Tl template. Functional images were rea- 
ligned, corrected for slice timing, coregistered with the 
structural volume, resampled to 2 ± 2 ± 2 mm, normal- 
ized into MNI space using parameters obtained from 
the structural normalization process, and finally 
smoothed with an 8 mm full width half-maximum 
Gaussian kernel to reduce spatial noise and to compen- 
sate for residual misregistration in the spatial normaliza- 
tion process. High-pass temporal filtering with a cutoff 
of 128 seconds was included in the first-level statistical 
model to remove the effects of low-frequency physiolo- 
gical noise. Serial correlations in the fMRI time series 
caused by aliased biorhythms were estimated using an 
autoregressive AR(1) model. 
First-level model 

After preprocessing, each subject's functional time series 
was modeled using a general linear model, with six 
regressors signifying the condition onsets and durations: 
the context, onsets of the CS+ and the CS-, the US, the 
omitted US, and the CS- offset. Movement parameters 
derived from the realignment step for x, y, z and roll, 
pitch, and yaw were included in the model to reduce 
the effects of residual motion-related noise. The experi- 
mental effects of interest were identified using a statisti- 
cal model containing boxcar functions representing each 
of the six experimental conditions, convolved with the 
SPM8 canonical hemodynamic response function. 

We focused our analysis on the responses to three of 
the experimental events: 1) the CS+, 2) the US and 3) 
the omitted US. The CS- cue and the interval following 
its offset were used as control contrasts. This design 
allowed us to investigate not only fear-conditioned 
responses but to also directly contrast brain responses 
induced by the actual delivery of the US with brain 
responses induced by the immediate expectancy and 
then the unexpected omission of the US. Because both 
the US and omitted US intervals occurred immediately 
after the offset of the CS+, a combination of anticipa- 
tory- and prediction-error components influencing the 
shock response could be identified in the omitted US 
contrast, as limited by the temporal resolution of the 
present fMRI paradigm. Multi-collinearity between 
regressors may lead to instability of parameter estima- 
tion and a consequent reduction in sensitivity. In the 



present design, the inherent multiple collinearity of the 
classic conditioning design (a cue followed by a shock) 
was not problematic for these three contrasts, as the 
correlations between the regressors were very low (r 2 < 
0.01). However, because of multiple collinearity between 
regressors identifying cue-induced activations and activ- 
ity immediately after the termination of the cue (r 2 = 
0.36), we did not examine this effect. 
Second-level model 

First-level contrast images representing the effects of CS 
+, CS-, US, omitted US and CS- offset were obtained for 
each subject, and modeled at the second level using a 
mixed-effects linear model with subject, group and task 
factors. At the second level, voxel-wise contrasts com- 
pared the PTSD group with the TENC group for CS+ 
versus implicit baseline, CS+ versus CS-, US versus 
implicit baseline, omitted US versus implicit baseline, 
and US versus omitted US. 

Recent studies of pain responses in PTSD have 
reported alterations in activity in the putamen, amyg- 
dala, hippocampus, cingulate, insula and caudate 
nucleus [11,12,23]. Based on these studies, we also per- 
formed restricted voxel-wise analyses within these 
regions defined from the Anatomical Automatic Label- 
ing (AAL) library [24,25]. Whole-brain and a priori 
regional activations of more than 10 contiguous voxels, 
adjusted for multiple comparisons within the search 
volume thresholded at P < 0.05 family-wise error 
(FEW), corrected for multiple comparisons, are 
reported. 

Symptom correlations 

First-level contrast images representing the effects of CS 
+ and US were obtained for all subjects with PTSD, and 
modeled at the second level using a linear regression 
model with CAPS scores. Correlated clusters of more 
than 10 contiguous voxels adjusted for multiple compar- 
isons at P < 0.05 FWE are reported. Correlation models 
such as these have been criticized for overinflated corre- 
lation coefficients [26], so we therefore did not plot cor- 
relations or report r values. 
Post hoc functional connectivity analysis 
To further understand alterations in US responses in 
PTSD, we performed a psychophysiological interaction 
(PPI) analysis, a technique allowing estimation of how 
inter-regional BOLD signal correlations vary with 
respect to an experimental psychological context. Based 
on the univariate results, we investigated whether the 
delivery of the US (versus omitted US) modulated the 
inter-regional functional connectivity between seed 
regions and targets elsewhere in the brain differently in 
PTSD and TENC subjects. Four seed regions (the left 
and right putamen, the left middle frontal gyrus, and 
the left parahippocampal gyrus) were chosen based on 
the peak observed hyperactivation to the US of these 



Linn man et al. Biology of Mood & Anxiety Disorders 201 1, 1:8 
http://www.biolmoodanxietydisord.eom/content/1 /I /8 



Page 5 of 1 2 



regions in the PTSD group. We contrasted inter-regio- 
nal influences of the seeds during the US interval with 
those during the omitted US interval to elucidate how 
the delivery of the US modulated seed functional con- 
nectivity to other brain regions. The corresponding PPI 
design matrix had three columns: the interaction 
between the experimental manipulation and seed time 
series, the main effect of experimental manipulation (US 
or omitted US), and the physiological effect (hemodyna- 
mically deconvolved seed region time series). In accor- 
dance with the task-related univariate analysis, we also 
included six estimated movement parameters to account 
for any residual effects due to inter-scan head motion. 
The voxel-wise regression of voxel time series from the 
entire brain on the three predictor variables resulted in 
a first-level PPI contrast for each seed region in every 
subject. Subsequently, these PPI contrast images were 
modeled in a second-level group analysis using a mixed- 
effects repeated-measures linear model with subjects, 
diagnosis and seed region as factors. This second-level 
analysis resulted in a statistical parametrical map of 
regions that displayed differences between TENC and 
PTSD subjects in alterations of seed functional connec- 
tivity depending on the presence or absence of the US. 
Clusters of more than 10 contiguous voxels adjusted for 
multiple comparisons at P < 0.05 FWE are reported. 

Results 

Shock levels and skin-conductance responses 

PTSD and TENC subjects chose similar shock levels and 
displayed similar levels of autonomic skin-conductance 
responses to the CS+ and the US. Subjects with PTSD 
displayed slightly higher skin-conductance responses to 
the non-reinforced CS- cue, but the magnitude of 
responses to the CS- was generally low (Table 2; see 
Additional file 1 Figure S3. Figure S4). There were no 
significant differences between PTSD and TENC sub- 
jects in brain activations to the CS+, the CS- or to the 
contrast CS+ versus CS-. 

Table 2 Shock levels and skin-conductance responses (in 
square root microSiemens, uS 1/2 ) to the conditioned 
stimulus, the unconditioned stimulus (US) and the 
omitted US ± standard error (significant P values in 
bold). 





PTSD a 




TENC b 




P 


Shock level, mA 


2.13 


± 


0.16 


2.20 ± 


0.18 


0.77 


CS C + response, |jS 1/2 


0.27 


± 


0.08 


0.17 ± 


0.04 


0.24 


CS- response, |jS 1/2 


0.08 


± 


0.04 


-0.01 ± 


: 0.04 


0.04 


Unconditioned response, |jS 1/2 


0.79 


± 


0.11 


0.71 ± 


0.07 


0.56 


Omitted US response, mS 1/2 


0.26 


± 


0.05 


0.21 ± 


0.03 


0.36 



a Post-traumatic stress disorder. 

b Trau ma -exposed normal control. 

C CS, conditional stimulus (positive or negative). 



Brain responses to the unconditioned stimuli 

For all subjects, the delivery of the US-elicited 
responses in the primary and secondary somatosensory 
cortices, insula, cingulate gyrus, thalamus and brain- 
stem (Figure 1). In the whole-brain between-group 
contrast, people with PTSD exhibited higher BOLD 
reactivity to the US in the bilateral putamen, left mid- 
dle frontal gyrus and left parahippocampal gyrus. 
Within the a priori regions of interest (ROI), subjects 
with PTSD displayed additional significant hyper-reac- 
tivity in the bilateral amygdala, left hippocampus, right 
dorsal ACC, bilateral posterior insula, left anterior 
insula, and bilateral caudate and putamen. No regions 
displayed lower reactivity in the PTSD group (Table 3; 
Figure 2). 

Brain responses to the omitted US 

During those trials in which the shock was unexpectedly 
omitted, at a time point when immediate anticipation, 
prediction error and relief from not getting shocked are 
intermingled in the hemodynamic response measurable 
in the current fMRI experiment, subjects with PTSD 
displayed significantly greater activity in hippocampus. 
This difference was driven by a lack of deactivation of 
the hippocampus to the omitted US in the PTSD group. 
No other significant between-group differences were 
noted in this analysis (Table 3). 

Brain responses to US versus the omitted US 

To partially separate sensory components of the US 
from anticipatory, prediction error and relief compo- 
nents, we contrasted US versus omitted US in PTSD 
versus TENC subjects (two-way ANOVA). Within the a 
priori ROI, this contrast showed significantly greater 
reactivity of the putamen and the anterior insula in the 
PTSD group (Table 3). 

Symptom correlations 

We found a significant correlation between PTSD symp- 
toms (CAPS scores) and US-induced BOLD activation 
in the right inferior frontal gyrus/anterior insula (MNIxyz 
= 28, 34, 0, cluster size = 329 voxels, Z = 7.93, P FWE cor _ 
rected = 0.002) (Figure 3). 

Functional connectivity of hyperactivated regions 

To further elucidate the findings from the univariate 
analysis, we performed post hoc functional connectivity 
analyses on the regions displaying higher BOLD 
response to the US in PTSD versus TENC subjects. We 
chose the US and the omitted US as the conditions pos- 
sibly inducing altered functional connectivity. In other 
words, we investigated how the functional connectivity 
of the seed regions was different during the delivery of 
the US compared with the omitted delivery of the US. 
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Figure 1 Main effect of unconditioned responses. Unconditioned stimulus (US) versus negative conditioned stimulus (CS)- offset in subjects 
with post-traumatic stress disorder (PTSD) and trauma-exposed non-PTSD control (TENC) subjects. Although the US was 'highly annoying' rather 
than painful, it elicited massive neuronal responses in classic pain regions including the brainstem, the thalamus, the contralateral sensory 
cortices, the insula and the middle cingulate. Activations are displayed at a threshold of T > 4 on a template magnetic resonance imaging (MRI) 
scan. 

L J 



Individual seeds of the putamen, middle frontal gyrus 
and parahippocampal gyrus were extracted and used in 
a psychophysiological interaction (PPI), that is, connec- 
tivity, analysis. 
Left putamen 

The left putamen displayed the most prominent hyper- 
activation in the PTSD group, consistent with previous 
studies using painful levels of heat [11,12], but its role 
in PTSD is less well described. We thus chose the puta- 
men region post hoc to perform a functional connectiv- 
ity analysis. The seed consisted of a 6 mm sphere 
centrered at the peak location (MNI xyz = -22, 0, -6). 
When comparing PPI effects of the left putamen 
between PTSD and TENC subjects, there were signifi- 
cant differences in functional connectivity shifts. The 
TENC group showed a a more positive PPI effect 
between the left putamen and the left temporal lobe 
including the amygdala. In other words, the shock deliv- 
ery led to a larger increase in functional connectivity 
between the putamen and the amygdala in the TENC 
group as compared to connectivity shifts in the PTSD 
group (Table 4; Figure 4a). 
Right putamen 

The right putamen seed (a 6 mm sphere) was centred at 
the peak location (MNI xyz = 18, 10, -12). Contrasting 
PPI effects of the right putamen between PTSD and 
TENC subjects showed differences in functional connec- 
tivity shifts similar to the left putamen seed, with a 
more positive PPI effect in the TENC group between 
the right putamen and the left temporal lobe including 
the amygdala, that is, indicating a larger increase in 
functional connectivity between these two regions at the 
shock delivery in the TENC group compared with the 
PTSD group (Table 4). 



Left middle frontal gyrus 

The left middle frontal gyrus seed (a 6 mm sphere) was 
centred at the peak location (MNI xyz = -22, -4, 48). 
Contrasting PPI effects of the left middle frontal gyrus 
between PTSD and TENC subjects showed no signifi- 
cant differences in PPI effects. 
Left parahippocampal gyrus 

The left parahippocampal gyrus seed (a 6 mm sphere) 
was centred at the peak location (MNI xyz = -18, -38, -4). 
Contrasting PPI effects of the left parahippocampal 
gyrus between PTSD and TENC subjects showed differ- 
ences in functional connectivity shifts, with a more posi- 
tive PPI effect in the TENC group between the left 
parahippocampal gyrus and right inferior parietal lobule, 
and the left medial frontal gyrus and the right precentral 
gyrus, that is, indicating a larger increase in functional 
connectivity between the seed and these three regions at 
the delivery of the shock in the TENC group compared 
with the PTSD group (Table 4 and Figure 4b). 

Discussion 

The predictable and 'annoying' shock (the uncondi- 
tioned stimuli (US)) elicited activation in several brain 
regions associated with pain processing [6], suggesting 
that stimuli may not need to be perceived as overtly 
painful to activate the so called 'pain matrix. Subjects 
with PTSD had significantly higher BOLD reactivity in 
the putamen, middle frontal gyrus and parahippocampal 
gyrus compared with the TENC subjects. The ROI ana- 
lysis further revealed the amygdala, hippocampus, dorsal 
ACC, insula and caudate nucleus to be hyperactive in 
the PTSD US response. These findings are consistent 
with previous studies indicating higher putamen and 
insula reactivity to heat pain in PTSD [11,12], and 
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Table 3 Univariate analysis results 



Contrast 


MNI b , peak x, y, z 


Z- sco re 


Size, voxels 


uusier r vaiue, rvvt 


Region 


ROI d 








US d 












PTSD e > TENC f 












Whole brain 


-22, 0, -6 


5.51 


1 1 1 


< 0.0001 


Left putamen 




(-24, 10,-14) 


-5.47 










(-32, 2, -4) 


-4.93 










18, 10-12 


5.47 


90 


< 0.0001 


Right putamen 




(20, 1 0, -4) 


-5.13 










-22, -4, 48 


5.25 


19 


0.003 


Left middle frontal gyrus 




-18, -38, -4 


5.07 


11 


0.007 


Left parahippocampal gyrus 


Putamen 


16, 12, -10 


5.34 


359 


< 0.0001 


Right ventral putamen 




(30, -20, 6) 


-3.79 










(28, -2, -6) 


-3.73 










-26, 2, -4 


5.1 1 


840 


< 0.0001 


Left putamen 




(-22, 8,-10) 


-4.88 










(-32, -12, -6) 


-4.75 








Amygdala 


-22, 6, -16 


4.3 


48 


0.01 


Left amygdala 




24, 2, -16 


4.27 


81 


0.005 


Right amygdala 


Hippocampus 


-18, -32, 4 


4.43 


36 


0.029 


Left hippocampus 




(-22, -24, -8) 


-3.74 










-26, -34, 8 


3.79 


18 


0.028 


Left hippocampus 




(-14, -32, 10) 


-3.52 










-26, -12, -12 


3.69 


25 


0.042 


Left hippocampus 




(-20, -4, 1 2) 


-3.62 








Cingulate 


8, 1 2, 42 


3.8 


19 


0.013 


Right dorsal acc 


Insula 


34, -20, 4 


4.97 


36 


0.006 


Right posterior insula 




-26, 14,-16 


5.05 


28 


0.008 


Left piriform insula 




-34, -6, 16 


4.14 


13 


0.017 


Left middle insula 




-36, -22, 4 


4.08 


17 


0.014 


Left posterior insula 


Caudate 


-14, 16, -8 


4.2 


67 


0.006 


Right caudate 




(16, 18, 0) 


-3.81 










16, 0, 16 


3.86 


24 


0.021 


Right caudate 




-12, 10, 14 


4.41 


297 


0.001 


Left caudate body 




(-18, -2, 22) 


-4.16 








TENC > PTSD 


No significant differences in whole brain or ROI 






Omitted US 












PTSD > TENC 












Whole brain 


No significant differences in whole brain 






Hippocampus 


22, -10, -20 


3.77 


54 


0.009 


Right anterior hippocampus 




(30, -10, -18) 


-3.46 








TENC > PTSD 


No significant differences in whole brain or ROI 






US minus omitted US 












PTSD > TENC 












Whole brain 


No significant differences in whole brain 






Putamen 


20, 12,-10 


3.89 


69 


0.003 


Right ventral putamen 




20, 0, 12 


3.48 


13 


0.032 




Insula 


-28, 16,-18 


4.05 


52 


0.016 


Left orbital insula 



Clusters are in bold, with sub-peaks in brackets. 

a Region of interest 

b Montreal Neurological Institute. 

c Family-wise error. 

d US, unconditioned stimulus 

e Post-traumatic stress disorder. 

f Trau ma -exposed normal control. 
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x = -18 



y = 12 



z = -8 




Shock response PTSD > TENC 



Figure 2 Increased putamen unconditioned stimulus (US) responses in subjects with post-traumatic stress disorder (PTSD). US 

responses in PTSD minus trauma-exposed non-PTSD control (TENC) displayed at T > 4 on a template magnetic resonance imaging (MRI) scan. 
The bilateral putamen displayed the most prominent hyperactivation in PTSD, no region displayed hypoactivation in PTSD. 



increased caudate signaling to heat pain after traumatic 
memory induction [23]. Subjects with PTSD displayed 
hyperactivation of the amygdala and dorsal ACC to the 
US, findings that are consistent with several previous 
studies suggesting amygdala and cingulate involvement 
in processing aversive stimuli [5,27-29] and hyperactiva- 
tion in PTSD [30]. Thus, a highly annoying but not 
overtly painful electrical shock can be used to replicate 
previous studies indicating alterations in PTSD pain 
processing. The benefits of using experimental stimuli 
that are not overtly painful are obvious. 




Figure 3 Post-traumatic stress disorder (PTSD) symptom 
severity correlates with the magnitude of unconditioned 
stimulus (US) reactivity. Positive correlations between blood 
oxygenation level-dependent (BOLD) US signal in the right anterior 
insula/inferior frontal gyrus and Clinician-Administered PTSD Scale 
(CAPS) score in the 19 subjects with PTSD investigated. The 
correlation map is displayed at a threshold of T > 4 on a template 
magnetic resonance imaging (MRI) scan. 



The CS+ led to comparable increases in skin conduc- 
tance in both the PTSD and the TENC group, suggest- 
ing no differences in the fear induced by the 
conditioned cues. Moreover, there were no differences 
in the neural response to the CS+, replicating behavioral 
studies that suggested intact fear learning in PTSD 
[31-35], but other studies have disagreed [36-38]. Auto- 
nomic responses to the US were of comparable magni- 
tude in PTSD and TENC subjects, as has been 
previously reported [34,36]. 

Subjects with PTSD displayed increased putamen 
reactivity to the US, an effect that was not observed 
after the omitted US. Moreover, the putamen hyperacti- 
vation to the US in the PTSD group remained signifi- 
cant when controlling for responses to the omitted US, 
(PTSD (US vs. omitted US) vs. TENC (US vs. omitted 
US); Table 3). We therefore interpret the putamen reac- 
tivity as driven mostly by sensation rather than altera- 
tions in anticipation of the aversive stimulus or a signal 
of error prediction. It has been shown that learning of 
the CS-US relationship allows healthy people to engage 
in endogenous pain inhibition [39-41]. When endogen- 
ous opioid neurotransmission is blocked with naloxone 
in healthy subjects, a moderately painful US leads to 
increased putamen and insular responses [42], similar to 
the present findings. There is evidence that people with 
PTSD have an altered opioidergic system [43-45], and 
the findings in the present study seem to support this. 

With regard to the functional connectivity results, two 
major effects were seen when comparing connectivity at 
the omitted US and the delivered US: TENC subjects 
(compared with subjects with PTSD), displayed a larger 
increase in putamen to temporal lobe/amygdala connec- 
tivity at the delivered US compared with the omitted 
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Table 4 Psychophysiological interaction analysis 



Seed 


MNI a , peakx, y, z 


Z-score 


Size, voxels 


Cluster P value, FWE b 


Region c 


Left putamen 












TENC d > PTSD e 


-44, -12, -16 


4.65 


295 


0.01 


Left temporal lobe 




(-32, 0, -30) 


4.01 






Left parahippocampal lobe/amygdala 




(-30, -10, -22) 


3.72 






Left parahippocampal gyrus 


PTSD > TENC 


No significant differences 








Right putamen 












TENC > PTSD 


-30, 0, -32 


4.54 


248 


0.019 


Left temporal lobe wm/amygdala 




(-48, 8, -30) 


3.97 






Left middle temporal gyrus 




(-56, 4, -26) 


3.44 






Left middle temporal gyrus 


PTSD > TENC 


No significant differences 








Left middle frontal gyrus 












TENC > PTSD 


No significant differences 








PTSD > TENC 


No significant differences 








Left parahippocampal gyrus 












TENC > PTSD 


44, -34, 42 


4.45 


449 


0.001 


Right inferior parietal lobule 




48, -38, 50 


-4.27 






Right inferior parietal lobule 




52, -26, 50 


-3.77 






Right postcentral gyrus 




-10, 66, 10 


4.42 


263 


0.017 


Left medial frontal gyrus 




(-14, 62, 0) 


-4.06 






Left medial frontal gyrus 




(-8, 62, 18) 


-3.45 






Left medial frontal gyrus 




34, -20, 68 


4.37 


487 


0.001 


Right precentral gyrus 




(34, -16, 46) 


-4.14 






Right precentral gyrus 




(22, -16, 74) 


-4.12 






Right precentral gyrus 


PTSD > TENC 


No significant differences 









a Montreal Neurological Institute. 
b Family-wise error. 

c Regions were identified based on nearest gray matter. 
d Post-traumatic stress disorder. 
e Trau ma-exposed normal control. 



US. Moreover, TENC subjects (compared with subjects 
with PTSD) displayed an increase in parahippocampal 
gyrus connectivity to the inferior parietal lobule, the 
precentral gyrus and the medial frontal gyrus (Figure 4a, 
b). We speculate that the PPI results may be indicative 



of an altered network processing of aversive stimuli 
PTSD, but further studies are needed to elucidate how 
the observed network connectivity changes might relate 
to the emotional and functional response to aversive sti- 
muli in PTSD. 



a 



x = -32 




x = -10 




Figure 4 Psychophysiological interaction (PPI) analysis. The PPI analysis contrasted the unconditioned stimulus (US) versus omitted US, 
modeling the change in regional coupling between the seed induced by the delivery of the US. The green sphere indicates the location of (A) 
the putamen seed and (B) the parahippocampal seed. Group differences (trauma-exposed non-PTSD control (TENC) > PTSD) in 
psychophysiological interaction effects are displayed at a threshold of T > 3 on a template magnetic resonance imaging (MRI) scan. 
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Higher right anterior insula US reactivity was corre- 
lated with higher PTSD symptom severity on the CAPS. 
This result is consistent with findings of Mickleborough 
and colleagues [23], who reported a positive correlation 
between CAPS scores and insula signaling during heat 
pain after traumatic memory activation. Dickie and col- 
leagues [46] reported a positive CAPS-insula correlation 
with remembered fearful faces. Felmingham and collea- 
gues [47] reported a positive correlation between CAPS 
and insula activation to masked fear faces. Thus, this 
relationship seems to be consistent whether the stimuli 
are perceived as painful, highly annoying but not pain- 
ful, encoded as aversive, or subliminally aversive. How- 
ever, Strigo and colleagues [12] found the rostral 
anterior insula reactivity to painful stimulation to be 
negatively correlated to CAPS avoidance scores. One 
possible explanation for this discrepancy is that in the 
Strigo experiment, pain stimulation was delivered with- 
out a predictive cue, whereas the study of Micklebor- 
ough and the present study entailed a predictive cue. 

The putamen is best known as a motor output 
region [48], but both nociception [49,50] and emotion 
[51] are processed in the putamen and the adjacent 
nucleus accumbens [52,53]. Moreover, the putamen is 
structurally and functionally connected to the medial 
orbitofrontal cortex and the amygdala [54], two regions 
implicated in PTSD. The observed increased putamen 
activation in PTSD corresponds to regions known to 
receive projections from the dorsal anterior cingulate 
and orbital cortex [55]. Altered basal ganglia function 
has previously been reported in PTSD for other para- 
digms: putamen blood flow and reactivity after script- 
driven imagery in PTSD has been correlated to aspects 
of flashback intensity and dissociative states [23,56], 
caudate and putamen volume are decreased in people 
with PTSD with headaches [57], and people with 
PTSD display diminished caudate and putamen signal 
to reward [58]. Thus, future studies could be aimed at 
elucidating PTSD alterations in both approach and 
avoidance behaviors involving corticostriatal circuitry 
[55,59]. 

Limitations and alternate interpretations 

We have interpreted the neuronal responses and con- 
nectivity modulations related to the omitted US as 
reflecting immediate expectancy (T am going to get 
shocked right now'). It could be argued, however, that 
these results also represent a prediction error and relief 
signal ('Oh, I guess I did not get shocked'). Our experi- 
mental design did not allow us to distinguish between 
these two responses, thus further studies with higher 
temporal resolution should be designed and conducted 
to further examine this point. 



Conclusions 

Previous studies have implicated the putamen and insula 
in altered pain processing in PTSD. The present study 
extends to this literature by showing that the functional 
activation of these brain regions is also altered in 
response to highly annoying but not overtly painful sti- 
muli, and that such alterations are likely to be driven by 
sensory aspects of the aversive stimulus rather than by 
expectancy- and prediction-error components. 

Additional material 



Additional file 1: Figure SI Paradigm design and timing. Figure S2 
Skin conductance responses within three different trials in one 
representative subject. Figure S3 Skin conductance responses to the 
conditioned stimulus. Figure S4 Skin conductance responses to the 
unconditioned stimulus. Figure S5 Skin conductance responses to the 
omitted unconditioned stimulus. 
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